An elastoplastic framework for granular materials 
becoming cohesive through mechanical densification. 
Part I - small strain formulation 

Andrea Piccolroaz, Davide Bigoni and Alessandro Gajo 

Dipartimento di Ingegneria Meccanica e 

Strutturale, Universita di Trento, 
Via Mesiano 77, 1-38050 Trento, Italia 
email: piccolroaz@ing.unitn.it, bigoni@ing.unitn.it, gajo@ing.unitn.it 

August 2, 2004 

Abstract 

Mechanical densification of granular bodies is a process in which a loose material 
becomes increasingly cohesive as the applied pressure increases. A constitutive 
description of this process faces the formidable problem that granular and dense 
materials have completely different mechanical behaviours (nonlinear elastic prop- 
erties, yield limit, plastic flow and hardening laws), which must both be, in a sense, 
included in the formulation. A treatment of this problem is provided here, so that 
a new phenomenological, elastoplastic constitutive model is formulated, calibrated 
by experimental data, implemented and tested, that is capable of describing the 
transition between granular and fully dense states of a given material. The for- 
mulation involves a novel use of elastoplastic coupling to describe the dependence 
of cohesion and elastic properties on the plastic strain. The treatment falls within 
small strain theory, which is thought to be appropriate in several situations; how- 
ever, a generalization of the model to large strain is provided in Part II of this 
paper. 



Keywords: Elastoplasticity; Granular materials; Mechanical densification; Forming; 
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1 Introduction 



1.1 A premise and the central problem 

Cold powder compaction is a process in which granular materials are made cohesive 
through mechanical densification. Subsequent sintering usually completes the treatment 
and yields the desired mechanical properties of the final piece. Since this process permits 
an efficient production of parts ranging widely in size and shape to close tolerances (Reed, 
1995), there is an evident related industrial interest. For instance, metallurgical (German, 
1984), pharmaceutical (Lordi and Cuitino, 1997) and forming of traditional (e.g. ceramic 
tiles and porcelain products) and structural (e.g. chip carriers, spark plugs, cutting tools) 
ceramics represent common applications. 

A crucial step in the above technique is the production of green bodies (i.e. the 
solids obtained after cold forming of ceramic powders) possessing enough cohesion to 
remain intact after mold ejection, being essentially free of macro defects, and handleable 
without failure in the subsequent treatments. This is not an easy task and indeed defects 
of various types are always present in the greens (Deis and Lannutti, 1998; Ewsuk, 1997; 
Hausner and Kumar-Mai, 1982; Thompson, 1981b), negatively affecting local shrinkage 
during sintering. In particular, defects can be caused by the densification process, which 
may involve highly inhomogeneous strain fields, or by mold ejection^. 

In the present article a new elastoplastic phenomenological model is formulated, de- 
veloped and implemented, which is capable of describing the deformation and consequent 
gain in cohesion of granular ceramic materials subject to mechanical loading. The model 
is calibrated by experimental data relative to a ready-to-press alumina powder, and nu- 
merical predictions are compared to experiments performed on simple forming processes. 
The proposed constitutive framework allows the possibility of simulating cold forming 
of pieces and to predict inhomogenities in cohesion, density, and elastic properties and 
residual stress distribution within the green, thus allowing a rational design of pieces. 

Models aimed at achieving results similar to those addressed here have been developed 
for metal powders (Ariffin et al. 1998; Brown and Abou-Chedid, 1994; Brown and Weber, 
1988; Khoei and Lewis, 1999; Lewis and Khoei, 1998; Lewis et al. 1993; Oliver et al. 
1996; Redanz, 1999; 2001; Redanz and Tvergaard, 2003; Sun and Kim, 1997) and for 
ceramic granular materials (Ahzi et al. 1993; Aydin et al. 1997a,b; Brandt and Nilsson, 
1998; 1999; Ewsuk et al. 2001; Keller et al. 1998; Piccolroaz et al. 2002; Zipse, 1997). 
Due to the plasticity of grains, metal powders are essentially different from ceramics, 
but nevertheless we point out that both the above classes of models cannot describe the 
features considered in the present papeiEl . In particular, the central point is that 

1 Glass and Ewsuk (1997) classify several types of damage such as for instance end and ring capping, 
laminations, shape distortions, surface defects, vertical cracks, and large pores. 

2 The model employed by Piccolroaz et al. (2002) is simply a modified Cam-clay calibrated on 
experimental results relative to ceramic powders. It is similar to many other available in the literature 
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the process of mechanical compaction requires the description of the transition 
from a granular to a dense or even a fully dense state. Since granular ma- 
terials are characterized by mechanical properties much different from those 
typical of dense solids, the constitutive modelling must describe a transition 
between two distinctly different states of a material. 

The above constitutive description has been obtained for the first time in the present 
article by suggesting a novel use of elastoplastic coupling theory (Hueckel, 1976), intro- 
ducing two micromechanically based hardening rules, and employing a recently proposed 
yield function (Bigoni and Piccolroaz, 2004). 

It may be also important to mention that, even if the emphasis is placed here on 
materials where cohesion increases with pressure, 'vice-versa', we believe our consti- 
tutive framework to be sufficiently general to describe progressive decohesion due to 
mechanically-induced damage of quasibrittle materials (concrete, cemented sand, and 
rock), a problem investigated by Lagioia and Nova (1995). 

Ceramic forming may involve deformations up to 50% and even greater, so that the 
need for a large strain formulation could be advocated. However, numerical simulations 
(some of which will be presented later) show that at least the 'gross' material behaviour 
is dominated by nonlinearities already occurring when deformations are small, so that 
the small strain formulation presented here should not be considered inappropriate. For 
completeness, nevertheless, a large strain formulation, requiring a new extension of the 
concept of elastoplastic coupling, is given in the Part II of this article (Piccolroaz et al. 
2005). 

1.2 The densification of ceramic powders and the related diffi- 
culties in the constitutive modelling 

The compaction of a ceramic powder is a process essentially consisting of three phases: 
(I) granule sliding and rearrangement, (II) granule deformation, and (III) granule den- 
sification. For obvious reasons, a sharp differentiation between the three phases cannot 
be established, but qualitative morphologies of granular arrangements during Phases I 
and II are shown respectively in the central and lower parts of Fig. [TJ with respect to the 
ceramic powder employed in our study (392 Martoxid KMS-96, defined in Appendix A 
and shown in the upper part of the figure in its loose state). 

The gain in cohesion starts with Phase II, which is marked by the so called 'breakpoint 
pressure', conventionally denoted by a change in inclination of the semi-logarithmic plot 
of density versus applied pressure (Fig. El referred to ceramic tablets obtained by uniaxial 
pressing of 392 Martoxid KMS-96 alumina powder). Since the densification process is 
often highly inhomogeneous, usually at least two phases coexist. Phase I always occurs 

and even if it can be considered perhaps sufficient for certain engineering purposes, it is far from realistic 
and does not address the major difficulties in modelling the densification process. 
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Figure 1: Micrographs of M KMS-96 alumina powder. The loose state is shown in the upper part, 
while granule arrangements corresponding to Phases I and II are shown in the central and lower part of 
the figure. 



in early volumetric deformation of granular materials, so that it has been thoroughly 
investigated for geomaterials and, as will be explained later, is characterized by an elastic- 
plastic behaviour in which the elastic response is nonlinearly dependent on mean stress. 
In contrast to Phase I, the other phases have been much less explored. The increase 
in cohesion^ becomes already substantial during Phase II, as shown in Fig. [3j where 

3 The degree of induced cohesion during mechanical densification depends on the applied pressure 
and is also strongly influenced by several factors, including stiffness, ductility and shape of the particles 
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Figure 2: Compaction diagram of M KMS-96 alumina powder. The breakpoint pressure (expressed in 
terms of axial stress during uniaxial strain), separating Phase I from Phase II, is denoted by p y . A few 
identical experiments have been performed, providing almost coincident curves, one of which is reported 
here. 



cohesion measured (on the same tablets employed for Fig. [2]) from biaxial flexure strength 
apparatus (defined by ASTM F 394) is reported versus the forming pressure (the axial 
compressive stress applied during uniaxial strain). 
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Figure 3: Biaxial flexure strength (following ASTM F 394) of green tablets from M KMS-96 alumina 
powder, as related to forming pressure under uniaxial strain in a cylindrical mold. Each point is the 
mean value taken over five tests, deviation was found to be negligible. 

The peculiar mechanism of variation in cohesion due to plastic deformation can be 
described by making recourse to the concept of hardening. In particular, we will assume 
that a yield function exists for a granular material, defining its elastic range, so that 
when the material is in the initial cohesionless state, the null stress state lies on the yield 

(Brown Abu Chedid, 1994). Experiments performed by D. Bigoni, G. Celotti, S. Guicciardi and A. 
Tampieri showed that penny-shaped particles produce a much higher increase in cohesion than spherical- 
shaped particles. 
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surface, Fig. HJ Now, if the material is subject to increasing hydrostatic compression, 
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Figure 4: Increase in cohesion as related to hardening. 



after an initial (small) deformation in the elastic range, an early development of plastic 
deformation occurs from a virgin state, corresponding to Phase I compaction. In this 
phase, the increase in cohesion is limited and almost negligible. However, when the 
pressure reaches the breakpoint value, so that material enters Phase II, the gain in 
cohesion becomes crucially important. In conclusion, to describe this process, we may 
employ a hardening law leading to a yield surface evolution of the type sketched in Fig. HJ 
where the yield surface shape distortion changes qualitatively, when the applied pressure 
p exceeds the breakpoint pressure p c b, expressed in terms of mean stress (with reversed 



At this level of description, one can get the impression that modelling the mechanism 
of cohesion increase during densification of granular materials could be pursued by simply 
employing an appropriate hardening rule. However, the elastic range of granular mate- 
rials cannot be properly described by linear elasticity so that the elastic response of the 
material changes, when the elastic range is modified during hardening. More in detail, 
during Phase I of densification, studies relative to geomaterials demonstrate definitively 
that the elastic law relating the volumetric deformation to the applied mean pressure is 
logarithmic 0, as sketched in Fig. |HJ Since the logarithmic law is simply not defined for a 
cohesionless material at null pressure, the increase in the cohesion of the material implies 
a modification to the elastic law, which becomes dependent on the plastic deformation, 
the physical quantity playing the role of the driving mechanism for densification. This 

4 The type of increase in cohesion could also qualitatively change after the Phase III of densification 
is entered, but since we do not possess enough experimental data relative to this behaviour (occurring 
however at very high pressures, not involved in the usual forming of ceramics), this is not accounted for 
in the modelling. We believe anyway that its consideration would be not difficult, once experimental 
results were made available. 

5 This law has reached the status of unchallenged law for fine grained soils, whereas it is strictly 
followed only during elastic unloading in coarse grained materials (Lambe and Whitman, 1969). 




6 



means that the elastic properties of the material must depend on the plastic deforma- 
tion, a feature that can be described by making recourse to the concept of elastoplastic 
coupling (Hueckel, 1976; Dougill, 1976). 




log (pressure) 



Figure 5: Logarithmic elastic law relating the elastic bulk modulus k to pressure (mean stress taken 
positive when compressive). 

In the following, we will introduce the three fundamental ingredients in the modelling 
of the densification processes, namely, (i) the yield function appropriate for the descrip- 
tion of the behaviour of granular materials, (ii) the nonlinear elastic model, coupled to 
plasticity, (hi) two micromechanically-based hardening laws. 

1.3 Summary 

Part I of this paper is organized as follows. After a brief explanation of notation, the 
constitutive laws are introduced in Sect. [2J together with the calibration to experimental 
results. Numerical simulations and comparisons to experimental results relative to simple 
forming processed are presented in Sect. [3j Extension of the constitutive model to large 
strain is deferred to Part II of this paper (Piccolroaz et al. 2005). 

1.4 Notation 

A standard, intrinsic notation is used throughout the paper, where vectors and second- 
order tensors are denoted by bold (the latter capital) letters. The scalar product, the 
trace operator and the transpose are denoted by the usual symbols, namely, 



AB = ti AB T , 




for every tensors A and B. The following invariants of stress er will be used 



P = 



tier 
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) 



(2) 
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where 6 G [0, 7r/3] is the Lode's invariant and 



J2 = -devcr • devcr, J 3 = -tr (devcr) 3 , devcr = cr —J, (3) 

in which devcr is the deviatoric stress and / is the identity tensor. 

We will employ two tensorial products between second-order tensors A and B, 
namely, 

(A®B)[C} = (B-C)A, (A®B)[C]= A ° + 2 ° B T } (4) 

so that I ®_I becomes the symmetrizing fourth-order tensor, defined for every tensor A 
as I~®I[A] = (A + A T )/2. 



2 The constitutive model 
2.1 The yield function 

During cold mechanical compaction of ceramic powders, the shape of the yield surface 
evolves from that typical of granular material (for instance, that of the modified Cam- 
Clay model) to that characteristic of a ductile, dense material (for instance, that of the 
Gurson model). To describe this process in terms of hardening, a suitably 'deformable' 
surface is needed. This was found by Bigoni and Piccolroaz (2004); in particular, the 
yield function takes the form 

q 

F((T,p c ,c) = f{p,p c ,C) + —r, (5) 

9(P) 

where p c and c are the parameters that will be assumed to depend on plastic deformation 
and thus they will define the hardening behaviour, q is the deviatoric invariant ([2"])?, f(p) 
is the function describing the dependence on the mean pressure p, eqn. assumed in 
the form 



-Mp cV /($ - $ m ) [2(1 - a)$ + a] if $ G [0, 1], 
f(p,Pc,c) = { ' (6) 

-00 i/$£[0,l], 

in which 

* = 4^, (7) 

Pc + C 

and g(8) describes dependence on the Lode's invariant 9 defined by eqn. §Z§ 3 taken as 

1 



cos[/3| — |cos _1 (7cos3# 



8 



The yield function described by eqn. flSJ-flEJ) has been motivated and explained in 
great detail by Bigoni and Piccolroaz (2004), where the interested reader is remanded 
for details. We mention here that the yield surface corresponding to eqns. (jHJ)- ([SJ) is 
extremely versatile and remains convex in a broad range of values of parameters M, p c , 

c, m, a, (3, 7. 

The yield function gradient Q = df/dcr, needed for the practical implementation of 
the model, is reported for completeness in Appendix B. 

2.1.1 Calibration of the parameter describing the yield function for alumina 
powder 

The problem now is to determine the yield surface parameters needed to describe the 
material considered here. For this purpose, we note that m and a, defining the shape 
of the meridian section, can be calibrated employing the experimental iso-density data 
in the (p, q) plane obtained from triaxial compression tests by Bonnefoy (as reported 
by Kim et al. 2002), under the assumption that iso-density curves correspond to yield 
surface sections (Kim et al. 2000; Kim et al. 2002). The values m = 2 and a = 0.1 were 
found to provide the best fitting, as shown in Fig. [6j where the data sets correspond to 
different levels of densification. 



300 




Figure 6: Meridian sections of the employed yield surface fitted to the iso-density data for alumina 
powder (taken from Kim et al. 2002) at different levels of densification (increase in cohesion is not 
appreciable at the scale of the figure). 

To our knowledge, experimental data are not available to define the deviatoric section 
of the yield surface for alumina powder (and more in general for ceramic powders). 
Therefore, parameters and 7 have been calibrated on the basis of the value of the angle 
of internal friction, determined for our alumina powder from a standard geotechnical 
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direct shear test apparatus^] as follows. 

The variation of the vertical displacement of the sample upper surface [a positive 
(negative) sign denotes a dilatant (contractant) behaviour] and of the applied shear force 
during shearing is shown in Fig. [7J 




Figure 7: Vertical (negative sign corresponds to contractant behaviour) vs. horizontal 
displacements (upper part) and shear force vs. horizontal displacement (lower part) of 
three samples of M KMS-96 alumina powder, for different vertical pressures (200, 500 
and 1000 kPa). 



It can be observed from the figure that the tested samples exhibited the typical 
behaviour of a loose sand, with compressive volumetric strains during shearing, without 
a softening phase. The fact that the samples sheared at 500 kPa and 1000 kPa of vertical 
pressure have the same volumetric strains is probably related to a slightly looser initial 

6 The apparatus consists of a shear box containing the sample, which is split in the middle height. 
When a normal force is applied, the horizontal force required to induce a movement of the upper half 
of the sample with respect to the lower half is measured. This test is useful for the evaluation of the 
friction angle of a granular material, like the alumina powder in Phase I of densihcation. The samples 
were formed by carefully pouring the ceramic powder within the shear box. Shearing was performed at 
a velocity of 0.2mm/min. 
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condition of the former sample. It can be observed that, except for the test performed 
under a vertical stress of 200 kPa, the steady state condition typical of the critical state 
is not reached and at the end of the test the strength and the volumetric strains of the 
samples are still slightly increasing. This effect is more pronounced at larger applied 
vertical pressures and is probably connected to the progressive deformation and rupture 
of the grains constituting the alumina powder occurring during shearing even at low 
confining pressures. This is consistent with the experimental evidences that very large 
shear strains are necessary to reach the steady state in sands, when grain crushing occurs. 

The maximum shear force reached at the end of the tests is plotted in Fig. [8] as a 
function of the applied vertical load. The results clearly lie on a straight line and may 
be interpreted following the Coulomb-Mohr failure criterion, to yield a friction angle 
approximately equal to 32°. Making recourse now to the Coulomb-Mohr model, the 
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Figure 8: Direct shear test on M KMS-96 alumina powder. 



ratio of the deviatoric section radius for triaxial extension, g(0), to that for triaxial 
compression, g(ir/3), is related to the angle of internal friction through 

fl(0) = 3 - sin 
g(n/3) 3 + sin0' U 

Eqn. ([9]) is not enough to determine the two unknowns (3 and 7. In the absence of further 
indications, parameter 7 was fixed to be equal to 0.9, providing a deviatoric section fairly 
close to the piecewise linear deviatoric section (corresponding to 7 = 1). From this value 
and employing eqn. (jUj), (3 = 0.19 follows. 

Parameter M, defining the pressure- sensitivity, has been calibrated making recourse 
to the concept of critical state, a peculiar state in which a granular material deforms at 
constant stress and constant volume. Since the critical state occurs at null plastic (and 
irreversible) volumetric strain rates, which (for the flow rule that will be specified later) 
is equivalent to trQ = 0, the expression for the yield function gradient eqn. (1B.ll) - (1B.3j) 
provides the following condition 

2(m + 1)(1 - a)$ m + ma$ m - 1 - 4(1 - a)$ - a = 0, (10) 
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which identifies the critical state. This implicit relation can be solved numerically, as 
soon as values for m and a have been selected, providing the value of $ corresponding 
to the critical state, 

= $*(m, a). 

In our case $* = 0.658 is obtained. The critical state point in the (p, q) plane is therefore 

P* = (Pc + c)$* - c, 

(11) 



q* = g(8)Mp cV /($* - $* m ) [2(1 - a)^ 



a\ 



Eqns. ffTTl) define, through parameter p c , the representation of the critical state line in 
the (p, q) plane. Due to the fact that c is a nonlinear function of p c , the critical state line 
results mildly curvecQ- Indeed, the line is straight only in the first phase of densification, 
as long as c = 0, and then deflects from linearity in the subsequent phase, when cohesion 
increases, approaching, after substantial plastic deformation and c ~ Cqo, again a straight 
line with the same slope as the initial line. In other words, the final inclination of the 
critical state line, reached after substantial increase in cohesion, is 

d S= VWE^WE^IR, (12) 

dp y 1 y 1 

It can be observed from eqn. ( TTTT) that, since c « p c throughout the densification 
process, the effect of c can be neglected, so that eqn. ( 1121) can be taken as the slope of 
the critical state line for the entire deformation process. With reference to the triaxial 
compression state, 9 = 7r/3, this slope is related to the angle of internal friction through 

dq _ 6sin0 
dp 3 — sin ' 

which gives a value of M equal to 1.1. 



2.2 Elastoplastic coupling 

Dependence of elastic response of a material on plastic deformation for describing degra- 
dation of elastic properties was suggested independently by Hueckel (1976) for soils and 
Dougill (1976) for concrete. The model was later developed by Hueckel and Maier (1977), 
Capurso (1979), Maier and Hueckel (1979), and Bigoni and Hueckel (1991). We will de- 
velop the concept of elastoplastic coupling in the way suggested by Bigoni (2000), which 
yields a symmetric constitutive operator in the specific case of associative flow rule. This 
fact follows from different constitutive assumptions, in particular, it will be assumed here 

7 Note that the definition of critical state is extended here to include cohesive and ductile materials; 
we assume for simplicity that cohesion does not change at critical state. 
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as in (Bigoni, 2000; Gajo et al. 2004) that the flow rule sets the so-called 'irreversible' 
strain rate, which in the case of coupling is different from the plastic strain rate. 

The necessity of elastoplastic coupling for modelling the densification process of gran- 
ular materials considered here may be motivated — as sketched in Fig. [9] — by the ob- 
servation that elastic unloading in a uniaxial deformation test shows a clear tendency 
toward a stiffening caused by the increase in cohesion. However, even if this effect would 




Displacement, mm 



Figure 9: Elastic stiffening during uniaxial deformation test (experimental results on M KMS-96 
alumina powder). 



be disregarded in a first approximation, elastoplastic coupling would always be needed, 
for the reason mentioned in the Introduction, namely, to match the increase in cohesion 
with the nonlinear, logarithmic elastic model usually accepted for granular materials. 

The basic concept of the elastoplastic coupling is that the elastic potential (p depends, 
in addition to the elastic, also on the plastic strain, so that 

--^ (e6 ' eP) (14) 



de e ' 

where e e and e p are the elastic and plastic components of deformation, respectively 
(obeying the usual addittive rule). 



2.3 Elastic potential for compaction during Phase I and early 
Phase II 

In the elastic range, granular cohesionless material obey the well-known logarithmic law 
sketched in Fig. |5l relating the increment in the void ratio Ae = e — (measured with 
respect to an initial value e ) to the current mean pressure 

Ae e = -Klog— , (15) 
Po 
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where the suffix e remarks that we are referring to the elastic range, po is the value 
of p corresponding to the initial void ratio eo, and k is the logarithmic bulk modulus. 
Eqn. (|T5|) is the starting point to obtain the nonlinear elastic potential employed in 
the Cam-clay model (Roscoe and Schofield, 1963; Roscoe and Burland, 1968), which 
is particularly suitable for the description of cohesionless granular media, or, in other 
word, during Phase I of compaction.]. However, our intention here is to describe the 
behaviour of materials which may increase (or decrease) cohesion as a function of the 
plastic deformation. Therefore, we have to introduce a modification in the elastic Cam- 
clay potential, first, to include a cohesion and, second, to make this dependent on plastic 
deformation. The easiest way to do this is to introduce a plastic-dependent cohesion in 
eqn. ([15]) . playing the role of a modification to the mean pressure 

Ae e = -/dog P + c( , 6P) . , (16) 

where c (e p ) is the (positive) parameter describing the cohesion and depending on plastic 
deformation. In particular, the cohesion is assumed to depend on the volumetric compo- 
nent only of plastic deformation. This may be motivated by qualitative micromechanical 
considerations. Following, for instance, the Rowe (1962) model of a granular material, 
a shear deformation yields a decrease (increase) of cohesion when accompanied by dila- 
tancy (contractivity), Fig. [101 Obviously, the relation ( Fl6|) is meaningful only in an early 




<j — dilatancy= loss of cohesion 

contraction = increase in cohesion 1> 

Figure 10: The mechanism of increase and loss of cohesion visualized in terms of the Rowe model. 



stage of Phase II compaction, since unrealistic small values of void ratio are predicted at 
increasing p. However, it is important to realize that the description of the initiation of 
Phase II is crucial, since the material goes from a cohesionless state to a solid state. 

8 The value of the logarithmic bulk modulus k, which governs the elastic behaviour of the material in 
the first phase of densification, was deduced from experimental results for the alumina powder considered 
here by measuring the slope of the curves obtained by loading and unloading in uniaxial strain tests. 
For this evaluation, we have assumed a constant ratio between the horizontal ah and vertical a v stresses 
equal to 0.47, as deduced from the formulae 

— = 1 — sm <p, 
a v 

which is usually employed for granular media (Jaky, 1944). 
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We are in a position now to proceed with eqn. f lT6|) in the way usually followed in 
the derivation of the Cam-clay model, i.e. in the case of eqn. (I15p. Therefore, assuming 
incompressibility of grains, the volumetric elastic deformation is given by (e — eo) / (1 + eo), 
so that eqn. ( Tl6|) defines a volumetric nonlinear elastic law, to be added to a linear elastic 
deviatoric constitutive equation. The resulting elastic potential is therefore 

(P(e e ,e p ) = -|(tre e ) 2 + ctre e + ^(p + c)exp (~^J +^e e -e e , (17) 

where /x is the elastic shear modulus and k — k/(1 + e ). 

The elastic potential ( TT7|) . with parameters fi and k taken independent of the plastic 
deformation, is certainly suitable to describe the behaviour of the material in Phase I of 
the compaction process, since here the material is still granular and c = 0, so that the 
model reduces to the celebrated Cam-clay model. Moreover, the elastic potential (JT7J) 
can describe the initiation of Phase II of densification, where the increase in cohesion 
starts to play a role. On the other hand, the potential ffTTj) embodies a nonlinear increase 
of elastic bulk modulus with mean pressure, a feature which is clearly unrealistic during 
the late Phase II of deformation. 



2.4 Elastic potential for compaction during the late Phase II 

In the late Phase II, the material becomes more similar to a porous solid than to a 
granular body so that its properties progressively change. Experimental evidence of the 
mechanical behaviour of green's bodies formed at pressures corresponding to Phase II 
deformation is scarce. In particular, experiments by Zeuch et al. (2001) on alumina 
powder evidence that elastic properties become linear functions of forming pressure p c . 
Therefore, an elastic potential such as 

4>&#) = (|-0(tre e ) 2 + Ai e e .e e , (18) 

can be expected, where both /i and K are linear functions of p c . Though eqn. (fl8l) 
is very simple, we believe that the introduction of more complicated laws would be 
straightforward. 



2.5 Elastic potential for compaction during Phases I and II 



During deformation, the material behaviour changes from that described by potential 
(TTTj) to that corresponding to potential (fl8|) . In order to describe this transition in the 
material behaviour, we modify the elastic potential (fl7|) as follows 



-(tre 



e\2 



ctr e 



(19) 



+ (Po + c) 



(tre 



e\2 



2k 



+ d 1/n kexp 



tr e e 



d}/ n i 
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where d > 1 is a 'transition' parameter depending on the plastic volumetric strain through 
the forming pressure and governing the passage from logarithmic to linear law of elastic 
bulk modulus. Note that for d = 1 and in the limit d — > oo the potentials (|17p and ( TTBj) 
are respectively recovered. Finally, n > 1 is a material constant defining the decay of the 
exponential term. Moreover, the elastic shear modulus /i is taken dependent on plastic 
volumetric strain through parameters d and c as follows 



/'('/) = /'(i + c ( d- 2 ) /'i- 



(20) 



where /io, and /ii are positive material constants. 

In conclusion, the nonlinear elastic stress/strain law may be obtained from eqn. (|T9|) 
and results dependent on the plastic strain through c and d (the dependence is often 
omitted in the following for conciseness) 



-/itr e e + c 



+ (po + c) 



1\ tre e 



— exp 



tr e e 



d 1 / n k 



Taking the time derivative^! of f[2"TT) . we get the rate equations 



I + 2fie e 
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(22) 



where c, d, and fi arise from elastoplastic coupling (c— d — fi — Oxn the usual uncoupled 
models) and the elastic fourth-order tensor E, together with its inverse (restricted 
to the space of all symmetric tensors) are given by 



E 



I ® I + 2/xJ® J, E 



-i 



2fj, - 3K t 



(23) 



in which the tangent bulk modulus K t depends on the plastic deformation through c and 
d and on the elastic deformation in the following way 



Po + c 



K 



d-- + d- 1/n exp 
d 



tr e t 



d 1 / n k 



(24) 



9 Time is intended here as a loading parameter, while the material behaviour is assumed inviscid. 
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We need now to specify the particular dependence of the transition parameter d on 
the forming pressure p c . To this purpose, we note that, assuming the existence of a 
saturation threshold for the value of the cohesion, the asymptotic behaviours of the 
bulk modulus K t , eqn. (121]) . and the shear modulus /i, eqn. (1201) . as d — > oo are 

K t ~ + oo ^ ^ ^ Coa^id, (25) 

respectively. Since experimental results by Zeuch et al. (2001) suggest that fi and !T t 
become linear functions of p c for large values of forming pressure, inspired by ( 1251) . we 
assume for simplicity that parameter d is a linear function of forming pressure p c , for 
values of pressure superior to the breakpoint threshold p c b, so that 

d = 1 + B < Pc - Pcb >, (26) 

where B is a positive material parameter and the symbol <> denotes the Macaulay 
brackets operator (defined for every scalar a as < a >= (a + \a\)/2). 

It has to be noted that the determination of parameters B, ^ and n is not easy. 
In principle, parameters B and \i\ can be obtained matching the asymptotic behaviour 
(1251) with experiments of the type performed by Zeuch et al. (2001). However, precise 
determination of elastic constants of green's bodies is certainly difficult, moreover, n 
needs also to be determined. We have estimated B, fix on the basis of the experiments 
by Zeuch et al. (2001), finding the values 5=0.18 MPa" 1 , ^=64. 

Regarding the constant n, we have plotted the evolution of the tangent elastic bulk 
modulus during a hypothetical isotropic compression test at different values of n. This 
has been possible employing eqn. ( 1241) with d given by eqn. ( 1201) and formula (1321) for the 
cohesion c that will be introduced later. In this way, K t depends on the forming pressure 
p c (which has been taken coincident with — trer/3) and on tre, which has been evaluated 
numerically solving the trace of eqn. (T2T1) . Results have been reported in Fig. [TTJfor n = 
{1, 6, 60, 600}, from which we note that for n — 1 unphysical behaviours appear, but these 
are not any more evident starting from n — 6. In the absence of experimental results, we 
have taken n = 6 in our analyses. We note also that sensitivity to this parameter becomes 
very low for values of n superior to 60, so that a refined determination of parameter n 
would require ad hoc experimental investigation. 

2.6 Two micromechanically-based hardening laws 

In order to further develop eqn. (T22]) . evolution laws for the hardening parameters repre- 
senting the forming pressure p c and the cohesion c are needed, providing the functional 
dependences of these parameters on the plastic deformation. At this point, recourse to 
micromechanical considerations and experimental evidence becomes necessary. 
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Fi gure 11: Simulated evolution of tangent elastic bulk modulus during isotropic compression. Four 
values of n are considered: n — 1, 6, 60, 600. 



2.6.1 The relation between the forming pressure p c and plastic volumetric 



Parameter p c is related to the plastic deformation and this relationship can be determined 
employing the micromechanical model proposed by Cooper and Eaton (1962). This takes 
into account the fundamental fact that compaction can be divided into the three phases 
mentioned in the introduction. Based on statistical micromechanics considerations and 
validated on several experimental results on ceramic powders, Cooper and Eaton (1962) 
provide a double-exponential law describing the first two phases of densification in terms 
of the relation between the plastic increment of void ratio Ae p and the pressure parameter 



where ai, 02, Ai and A 2 are material (positive) constants. In particular, coefficients —e^ai 
and — eoa2 denote the increment of void ratio that would be achieved at infinite pressure 
by each of the two processes of densification, so that < a\ + a-i < 1. Coefficients Ai 
and A2, having the dimension of stress, indicate the magnitude of the pressure at which 
the particular process of deformation has the maximum probability density. 

Assuming incompressibility of the grain constituents^, the plastic volumetric defor- 

10 Ceramic powders are usually obtained through spray-drying and are formed by granules with 
dimensions ranging between 50 and 200 /im, coated with the binder system. Granules are itself aggregates 
of crystals having dimensions on the order of 1 fim. The crystals are here assumed to be incompressible, 
while compressibility of granules is due to their internal voids. 



strain 




(27) 
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mation is related to the plastic void ratio increment according to the rule 



Ae p = (l + e )tre p , (28) 

we get from eqn. (J27j) 

tr e p = —a\ exp f -J — a 2 exp ^ -J , (29) 

where a« = e^il (1 + eo), i = 1, 2. 

The hardening rule (1291) is calibrated to describe uniaxial strain experiments, in which 
the permanent volumetric deformation has been measured at various forming pressure. 
In particular, taking the values Ai = 1.8 MPa, A 2 = 40 MPa, a\ = 0.37, a 2 = 0.12 
gives the excellent interpolation presented in Fig. [121 where the volumetric plastic strain 
tre p is reported versus the hardening parameter p c . It may be noted that we restrict the 
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Figure 12: Plastic volumetric strain as a function of the forming pressure p c . Experimental data 
relative to M KMS-96 alumina powder are compared with the model by Cooper and Eaton (1962) (solid 
line). 



attention to the double-exponential law (12"9"|) for simplicity, but it would be certainly not 
difficult to include more complicated relationships, which — as suggested by Cooper and 
Eaton — could include an arbitrary number of exponentials and therefore describe also 
the third phase of compaction behaviour. 

Eqn. (I29p defines an implicit relation between the plastic deformation tre p and the 
forming pressure p c , which becomes explicit in terms of rates. In particular, the rate of 
eqn. ( 129]) gives 

Pc = , ^ - n ~r e p , (30) 

OiAj exp I J + a 2 A 2 exp I 

V PcJ V Pc 
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providing the first hardening rule. 



2.6.2 Modelling the increase of cohesion 

Concerning the dependence of cohesion c on the forming pressure p c , we could recourse to 
models of adhesion between particles. We note that the celebrated JKR model (Johnson 
et al. 1971) (and also variants like for instance the DMT model) is not applicable in our 
case, since the ceramic granules considered here are highly plastic and the adhesion force 
cannot be treated as independent of the granule deformation. It seems more appropriate 
to recourse to the Bowden and Tabor (1950) model for adhesion between surfaces in 
contact. In this model, the cohesion depends on the real contact area and therefore on 
the normal pressure. Assuming a Herzian contact between spheres and that the cohesion 
c is a linear function of the contact area, the following relationship 

c~V 2l \ (31) 

is found, in which p is the contact pressure. The major concern with condition ( f3~li) is 
that it does not predict a limit for the increase of adhesion with pressure, which is a 
clear experimental evidence. More in detail, our experimental results reported in Fig. [3] 
can be manipulated employing our model and expressed in terms of cohesion c versus 
the forming pressure p c . The results are shown in Fig. [121 where interpolation using the 
Bowden and Tabor approximation (f3"T|) is also reported (dashed). The solid line in the 
figure is obtained employing the following law 

c = Coo [1 - exp (-r < p c - p cb >)] , (32) 

where p cb is the breakpoint pressure, Coo and V are two positive material parameters, the 
former defining the limit value of cohesion reached after substantial plastic deformation, 
the latter related to the 'velocity of growth' of cohesion. The values V = 0.026 MPa -1 , 
Cqo = 2.3 MPa, and p cb = 3.2 MPa have been found to provide an excellent interpolation 
to experimental data. Note that the Macaulay bracket is needed since the cohesion is 
null before the threshold value defined by the breakpoint pressure is attained. 
Taking now the rate of eqn. fl3"2"l) and using fl3"Uj) yields 



where 



c = £ 2 tre p , (33) 
c^T H(p c -p cb )exp[-T(p c -p cb )]p 2 c 

~ ; — ma — ; — ma' ( } 

eiiAi exp + a 2 A 2 exp 

V PcJ V PcJ 

providing the second hardening rule. Note that symbol H in eqn. (|34|) denotes the 
Heaviside step function (defined for every scalar a as H(a) = 1, if a > 0, H(a) = 
otherwise) . 
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Figure 13: Variation of the cohesion c as a function of the forming pressure p c . Experimental results 
relative to M KMS-96 alumina powder are compared to our model (|32|) (solid line) and to that by 
Bowden and Tabor (1950) (dashed line). 



2.7 The elastoplastic coupling in rate form 

Taking the time derivative of eqns. ( )26|) and (120]) and considering eqn. (130]) we get the 
dependence of elastic parameters on plastic strain in rate form 

<i = £ 3 tre p , and /i = ^4tre p , (35) 

where 

t _ BH( Pc - Pcb )p 2 c 



. Ai\ ( A 2 

aiAiexp ( -— ) + a 2 A 2 exp ( — — - 



Pc/ \ P, 

^={d~l) ^ + c ( 1 + ^) ^3- 

Employing the evolution equations (133]) and (135]) into the stress rate equations f]22]) . 
we may write 

& = E[e e ] + P[e p ], (37) 

where the fourth-order tensor P defines the contribution of the elastoplastic coupling, in 
the sense that P is null in the usual, uncoupled plasticity and is defined as 

P = ® / + 2£ 4 e e ® I, (38) 
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in which 



£ 5 = -g£4tre e + £ 2 



1\ tre e 



d 



K 



— exp 



tr e e 



K 



— - 1 

+ d 2 nd^+ l l' 



cxp 



d 1 / n n / 
tr e e 
d l l n k 



(39) 



when elastoplastic coupling occurs. 

Introducing now the 'irreversible' (in an infinitesimal stress cycle) strain rate 



in which 
and 



6 



e* = G[e p ], 
I® J - E _1 P = I® J + £ 6 I ® I + &e 



k 



3K t 9/iK t 

we may transform the rate equation (1371) into the equivalent form 

a = E[e] — E[e*]. 
Since tensor G is assumed positive definite, implying that 



e p > 0, 



the inverse of G is given by 



-1 _ 



where 



?8 



6 



J, 

6 



(40) 

(41) 
(42) 

(43) 

(44) 

(45) 
(46) 



1 + 3£ 6 + £ 7 tr e e 1 + 3£ 6 + £ 7 tr e e 

The irreversible deformation rate e l defined by eqn. (1401) is the rate of deformation 
which is not recovered in an infinitesimal stress cycle. This should not be confused 
with the plastic deformation rate e p , which can be only detected upon unloading at zero 
stress. This concept is illustrated in Fig. [J4] with reference to a hypothetical volumetric 
stress/strain law, where the rate deformation at unloading E _1 [<r] is also indicated. Due 
to the plastic increase in elastic stiffness, we note that (122]) implies that the plastic e p 
and inelastic e l rates are different. 
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2.8 Flow rule 

A crucial point is now the definition of the flow rule, that following Bigoni (2000) is 
assumed to govern the irreversible strain rate (instead of the plastic, see Gajo et al. 2004 
for a discussion on this issue) 

e l = XP. (47) 

Experimental evidence for granular material supports the use of a deviatoric associative 
flow rule, which is adopted here. Therefore, tensor P is related to the yield function 
gradient Q through 

p = q _ e(l- 3 *)t,O f[ 0£f£1 _ m 

where e is a parameter governing the entity of volumetric nonassociativity, so that e = 
gives the associative flow rule. All indirect evidences point to flow rule nonassociativity 
for ceramic powders, so that we feel that eqn. ( 148]) is appropriate; however, experimen- 
tal evidence for alumina powder is not available and associativity will be assumed for 
simplicity in the following numerical simulation. 

The rate constitutive equations can now be obtained via Prager's consistency, so 
that F = during plastic deformation. Imposing this condition suggests the following 
definition of hardening modulus 

1 fdF OF \ 
h = 7-1 — p c + —c ), (49) 



A \dp c dc 

which is positive in the case of hardening, negative for softening and null for ideally 
plastic behaviour. The derivatives of F with respect to the hardening parameters p c and 
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c appearing in eqn. (j4"9~j) are given by 



ur i 

— = -My/($ - $ m ) [2(1 - a)$ + a] 



(50) 



'c 



p c (p + c) (1 - m®" 1 - 1 ) [2(1 - a)$ + a] + 2(1 -«)($- $ m ) 



(Pc + c) 2 2^(§-<& m ) [2(l-a)$ + a] 



and 



<9F 

9c 



— M 



Pc(Pc - g) 

(Pc + cy ' 

(1 - m^™- 1 ) [2(1 - a)$ + a] + 2(1 - a) ($ - $ m ) 



(51) 



2y / ($-$ m ) [2(1 - a)$ + a] 



2.9 The final rate equations 

Employing definition (j4"9j) into Prager's consistency yields the elastoplastic rate equations 



It may be noted that the elastoplastic tangent operator becomes symmetric in the specific 
case of the associative flow rule, P = Q. 

3 Numerical simulations 

The proposed constitutive model was implemented into UMAT, the subroutine avail- 
able within the commercial finite element code ABAQUS (Ver. 6.3; Hibbitt, Karlsson 
& Sorensen, 2002, Pawtucket, RI, USA). The employed numerical integration scheme 
was the so-called 'cutting-plane algorithm', proposed by Simo and Ortiz (1985), Ortiz 
and Simo (1986), Simo and Huges (1987). A full Newton- Rapson scheme has been em- 
ployed for the solution of the nonlinear finite element problem. Parameters of the models 
employed for the simulations are summarized in Table [TJ 

3.1 Forming of tablets 

Simulations of uniaxial deformation of cylindrical samples, or, in other words, forming 
of tablets, are reported in Fig. [151 together with our experimental results, marked by 
spots. In the figure the applied vertical load (in kN) is reported versus the vertical 
displacement (in mm) and different forming pressures, equal to {60, 80, 100, 120} MPa, 
have been considered. The fact that the experiments are correctly simulated employing 




if F(er,p c ,c) = 0, 
if F((T,p c ,c) < 0. 



(52) 
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Table 1: Values of material parameters estimated from experiments for alumina powder 
392 Martoxid KMS-96. 



Material parameters defining the yield surface 
m = 2, a = 0.1, = 0.19, 7 = 0.9, M = 1.1 



Elastic logarithmic bulk modulus K = 0.04 



Material parameters defining the hardening rule (|29[) 
Aj = 1.8 MPa, A 2 = 40 MPa, 61 = 0.37, h 2 = 0.12 



Material parameters defining the hardening rule (|32 
T = 0.026 MPa" 1 , Coo = 2.3 MPa, p cb = 3.2 MPa 



Material parameters defining coupling rules (|2"6")l and (|2"0")) 



£ = 0.18 MPa~ 



6, = 1 MPa, /xi = 64 



Material parameter defining the flow rule e = 
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Figure 15: Experimental (data relative to M KMS-96 alumina powder) and simulated load vs. dis- 
placement curves for tablets formed at various final pressures: 60 (a), 80 (b), 100 (c), and 120 MPa (d). 
Different values of n have been considered. 



the finite element discretisation should be considered as a succesful feedback on the 
'consistency' of the model and its implementation. Moreover, we note that the model 
describes the progressive increase in elastic stiffness (visible at unloading) in qualitative 
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and quantitative agreement with experimental data. 
3.2 Forming of a simple ceramic piece 

Numerical simulations were performed to describe forming of the (axisymmetric) piece 
geometrically described in Fig. [16] and shown as a green body after forming in Figs. [16j 
In particular, four pieces have been formed at a final mean pressure of 100 MPa, 
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Figure 16: Geometry of the formed piece (dimensions in mm). 



Figure 17: Photograph of the formed green piece (5 g of M KMS-96 alumina powder has been used 
and a final mean pressure of 100 MPa has been reached). 

starting from 5 g of the M KMS-96 alumina powder employed in all our experiments. 
The axisymmetric mesh, employing 4-node elements (CAX4), used in the simulations is 
shown in Fig. [THJ 

The following assumptions have been introduced to simulate the entire forming pro- 
cess: 

• the die is rigid; 
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Figure 18: Initial mesh. 



• friction is neglected at the contact between powder and die walls; 

• the initial configuration is that shown meshed in Fig. [T8j 

After the initial state, defined by given initial values of void ratio and confining pres- 
sure (e =2.129 and p =0-063 MPa have been assumed, respectively), has been prescribed, 
the loading history is assigned, in terms of the following three sequential steps: 

1. forming is prescribed by imposing the motion of the upper part of the boundary 
(3.78 mm, corresponding to the value measured during forming at the final load of 
50 kN); 

2. unloading is simulated by prescribing null forces on the upper part of the boundary; 

3. ejection is simulated by prescribing null forces on all the boundary. 

The final deformed mesh (at the end of step 3), is reported in Fig. [191 superimposed 
on the initial mesh. It can be noted from the figure that the elements near the corner 
of the punch are excessively distorted so that results in this zone should be considered 
unrealistic. Comparing the meshes, it can be observed that the deformation suffered by 
the piece is moderately large. 

The hydrostatic stress component p (taken positive when compressive, upper part in 
the figure), the Mises stress q (central part in the figure) and the void ratio e (lower part 
in the figure) are reported in Fig. l2"Ulat the end of step 1, in Fig. EH at the end of step 2, 
and in Fig. [22] at the end of step 3. 

Regarding Fig. [20] (representative of step 1) we may note that, excluding the small, 
unrepresentative zone near the corner of the punch, the hydrostatic stress p ranges be- 
tween 40 and 100 MPa and the Mises stress q between 25 and 120 MPa, evidencing a 
high stress inhomogeneity. 
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Figure 19: Initial and deformed (step 3) meshes. 



Considering Fig. |5T] (representative of step 2) it may be important to note that 
residual stress is definitely high, due to the lateral constraint (provided by the walls 
of the mold) still present at the end of step 2. The knowledge of the lateral stress is 
important for practical purposes since the force needed for the ejection of the final piece 
can be estimated from this value employing the Coulomb friction law. In particular, a 
rough, but simple evaluation can be immediately obtained from numerical output at the 
end of step 2 employing the formula 

ejection force = atan0 (mean lateral stress x lateral surface of the piece), 

where is the powder friction angle (equal to 32° in our case) and a is a coefficient 
dependent on the roughness of the die wall and ranging between and 1. 

Regarding Fig. [22] (representative of step 3), we note that p and q represent, in terms 
of hydrostatic stress and Mises stress components, the residual stress distribution in the 
green body at the end of forming. 

Excluding again the unrepresentative zone near the corner of the punch, the hydro- 
static stress p ranges now between -1.5 and 6 MPa and the Mises stress q between 1 and 
6 MPa. Moreover, the void ratio varies between 0.6 and 0.9. It can be noted that the 
minimum void ratio is not associated with the maximum residual mean stress, it is rather 
associated with the maximum mean stress reached during loading (step 1). The results 
suggest also that two oblique zones of material are formed, the outer of which is subject to 
high compressive mean stresses, whereas the inner is subject to tensile stresses, creating 
a sort of truss-like mechanism. This can be considered representative of a situation where 
the tensile stresses tend to open possible microcracks, leading to defects formation in the 
green. It is however worth remembering that, even in the absence of macro defects, the 
mechanical behaviour of the green and the shrinkage during future sintering are deeply 
affected by the inhomogeneities in the residual stress and density distributions. 
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Figure 20: Hydrostatic stress p (MPa, upper part), Mises stress q (MPa, central part), and Void ratio 
e (lower part) distributions at the end of step 1 (after loading). 
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Figure 21: Hydrostatic stress p (MPa, upper part), Mises stress q (MPa, central part), and Void ratio 
e (lower part) distributions at the end of step 2 (after punch removal). 
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Figure 22: Hydrostatic stress p (MPa, upper part), Mises stress q (MPa, central part), and Void ratio 
e (lower part) distributions at the end of step 3 (in the green piece). 
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From the comparison reported in Fig. [19] between the initial mesh and that at the end 
of step 3, it becomes now possible to evaluate the springback effect. In contrast to the 
prediction of the simple model employed by Piccolroaz et al. (2002), it can be noted that 
now the model correctly predicts that the springback effect and the shape distortion are 
very small. In particular, the final diameter of the piece is estimated to be 0.1 mm larger 
than the inner diameter of the die, in agreement with our experimental observations. 

The cohesion c attained by the material at the end of the overall process is shown 
in Fig. 1231 upper part, whereas the elastic properties of the final piece are reported in 
Fig. [23] (central and lower parts, respectively), in terms of tangent bulk modulus K t and 
shear modulus fi. The inhomogeneity of the elastic properties evidenced in Fig. [23] is 
obviously a consequence of elastoplastic coupling. This effect and also the increase in 
cohesion have been not modelled in the simple analysis presented by Piccolroaz et al. 
(2002). It may be observed by comparing the maps shown in Fig. [23] with the map of 
hydrostatic stress p at the end of step 1 (Fig. [20] upper part) that there is a strong 
relation between mechanical properties gained by the material in the final piece and the 
mean stress reached during loading. This results from our analyses to represent the most 
important parameter in the entire forming process. 

Experimental and simulated load displacement curves during forming (natural and 
semilogarithmic representations are reported) are compared in Fig. [2H together with the 
results obtained by Piccolroaz et al. (2002) included (dashed) in the figure. Beyond the 
excellent agreement, it may be instructive to compare the present numerical simulation 
with that performed by Piccolroaz et al. (2002) employing a simple model. First, we note 
that the simple model approximately describes two straight lines, whereas the present 
model describes a curved line, much closer to the experimental results. Second, we may 
speculate on the limits of the approach presented by Piccolroaz et al. (2002) and of many 
similar models available in the literature; in particular, the modelling can be accurate 
enough, if intended to predict a 'global' force-displacement curve like that reported in 
Fig. [23] On the other hand, the simple model does not predict increase in cohesion and 
dependence of elastic properties on plastic deformation, so that the internal stress and 
strain distributions result almost completely different. 

A photograph of the lower side of the formed pieces is shown in Fig. (25] where we 
can note the formation of annular dark zones, evidencing a sliding between material and 
mold. This sliding is indeed predicted by the simulations, so that the radial displacement 
at the mold contact is reported in Fig. [2B1 superimposed to the photograph of one of 
the pieces shown in Fig. [25] so that, since the dark zone corresponds to the peak of the 
radial displacement, simulations again agree with experimental observations. 

3.3 The effects of large strains 

In order to have an insight on the possible effects of large strains, a few analyses have 
been performed using the option NLGEOM available on ABAQUS, still employing the 
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Figure 23: Distribution of cohesion c (MPa, upper part), bulk modulus K t (MPa, central part), and 
shear modulus fj, (MPa, lower part) at the end of step 3 (in the green piece). 
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Figure 24: Experimental (M KMS-96 alumina powder) and simulated load vs. displacement curves 
during forming of the piece shown in Fig. 1171 in natural and semilog representations. Results by 
Piccolroaz et al. (2002) are also reported (dashed). 



presented small-strain formulation. A representative result is shown dashed in Fig. |271 
relative to the forming of a tablet [see Fig. [15] (d)]. Obviously, a rigorous analysis 
would require a complete model reformulation and reinterpretation of experimental results. 
However, we believe that due to the fact that the deformations are only moderately 
large and do not involve rotations, the simulation shown in Fig. |27] should at least give 
some understanding of the main differencies between predictions expected from the two 
formulations. We note that the effect of geometrical nonlinearities yields, as expected, 
a stiffening of the response, but does not change the results qualitatively. However, the 
quantitative difference may be enough to suggest the interest in a large strain formulation, 
which is given in the Part II of this paper (Piccolroaz et al. 2005). 

4 Conclusions 

Results presented of this paper provide a constitutive framework to realistically describe 
forming processes of ceramic materials. Even if the experimental results are still in- 
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Figure 25: Photograph of the lower side of the formed pieces. 




Figure 26: Simulated radial displacements superimposed to the photograph of the lower surface of the 
formed piece. 

complete and the employed elastoplastic model has been developed in a small strain 
formulation, it has been shown that it is possible to predict: 

• the force needed for mold ejection, 

• the springback effect and related shape distortion of formed pieces, 

• the residual stress distribution, 
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Figure 27: Experimental (M KMS-96 alumina powder) and simulated load vs. displacement curves for 
uniaxial strain. The prediction obtained employing the option NLGEOM has been reported (dashed) 
together with the prediction relative to the small strain assumption. 



• the gain in cohesion and the final elastic properties, 

• the density distribution and the related presence of defects in the green body. 

The last of the above points is related to the prediction of defects in the sintered piece 
and therefore its investigation has an important consequences in the design of ceramic 
elements. 

In closure, we mention that the present modelling can be extended in different direc- 
tions. Introducing thermoplastic effects, the sintering phase might be covered by mod- 
elling, so that simulation could be extended to the entire production process. Moreover, 
both sintering aids and powder characteristics might enter the elastic-plastic constitutive 
laws, so that the optimal powder composition and morphology could be predicted for 
different forming problems. 
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APPENDIX A. Powder characteristics 



Calibration of the model has been performed on the basis of experiments both already 
available and carried out on a commercial ready-to-press alumina powder (96% purity), 
manufactured by Martinswerk GmbH (Bergheim, Germany) and identified as 392 Mar- 
toxid KMS-96. The data presented by the manufacturer are given in Table [2 It can be 
noted from the upper part of Fig. [I] that the granules have a mean diameter of 250 /jm. 



Table 2: Granulometric and density properties of the tested alumina powder. 



MWM 28 Vibration sieving 
sieve residue > 300 /im 
sieve residue > 150 /im 
sieve residue < 63 /im 


3.9% 
56.3% 
2.5% 


Bulk density (g/cm 3 ) 


1.219 


Green density (p = 50MPa) (g/cm 3 ) 


2.39 


Fired density (T=1600°C, 2h) (g/cm 3 ) 


3.77 
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APPENDIX B. Yield function gradient 



The gradient Q of the yield function (jHJ) is [a detailed derivation can be found 
(Bigoni and Piccolroaz, 2004)] 
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